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We report a ^''Cu NMR investigation in detwinned YBa2Cu30j, single crystals, focusing on the 
highly underdoped regime {y = 6.35 — 6.6). Measurements of both the spectra and the spin-lattice 
relaxation rates of ^"^Cu uncover the emergence of static order at a well-defined onset temperature 
To without a known order parameter as yet. While To is rapidly suppressed with increasing hole 
doping concentration p, the spin pseudogap was identified only near and above the doping content 
at which To — >■ 0. Our data indicate the presence of a critical hole doping pc ~ 0.1, which may 
control both the static order at p < pc and the spin pseudogap a.t p > pc- 



The superconducting coppcr-oxidcs (cupratcs) in the 
underdoped regime feature unusual states of matter, such 
as a pseudogap (PG), density- wave order (stripes), and 
the coexistence of magnetism and superconductivity. De- 
bates regarding the origin and the precise nature of those 
phases or related phenomena like a reconstruction of the 
Fermi surface at a quantum critical point are still ongoing 
actively [iHl]. In particular, interest in the underdoped 
YBa2Cu30j, (YBCOy) has been revived in recent years, 
during which significant progress in understanding those 
subjects has been made through experimental observa- 
tions of an electric liquid crystal (ELC) or nematic phase 
0, Q , quantum oscillations above a critical doping p-Q , 
and field- induced charge stripe order [l^l . 

Motivated by the recent literature and by the avail- 
able high-quality single crystals of underdoped YBCOy, 
we carried out a ^-^Cu NMR study of YBCOy to elucidate 
the underlying physics in the highly underdoped region 
of the compound on a microscopic level. While nuclear 
magnetic resonance (NMR) is a powerful local probe, so 
far, the majority of the NMR studies on the planar Cu in 
YBCOy has been performed on nearly optimal or slightly 



underdoped regions [ll|, ll2| , largely due to strong mag- 



netism which causes complicated static and dynamic ef- 
fects on NMR parameters. In this Letter, we show that 
a critical hole doping pc exists in the p-T phase diagram 
of YBCOj, beneath the superconducting (SC) dome at 
p ~ 0.1, below which a static order sets in and above 
which a spin pseudogap (PG) opens up in the low-energy 
spin excitation spectrum. 

The growth and characterization of detwinned YBCOy 
single crystals are described in Rcfs. [l3,[l3|- The single 
crystals investigated here have y = 6.35 (Tc = 10 K, 
p = 0.062), 6.4 {Tc = 21 K, p = 0.075), 6.45 (T^ = 35 K, 
p = 0.082), 6.5 (sample 1 with short ortho II correlation 
length [SJ]: Tc = 53 K, p = 0.106; sample 2 with long 
correlation length (~ 100 A): = 61 K, p = 0.114), 
and 6.6 {Tc = 61 K, p = 0.135), where p is the hole 
concentration per planar Cu determined from the c-axis 
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FIG. 1: ^^Cu NMR spectra in YBCOy at 7 T along the c axis 
for y = 6.35 (left panel) and 6.45 (right panel). Cu(l) denotes 
the Cu sites in the CuO chains, where the subscripts and 
2 represent the number of neighboring oxygen ions along the 
chain, and Cu(2) (dotted lines) denotes the planar Cu sites 
[Hi. An emerging new Cu spectral feature spO consisting of 
two resonance lines indicated by two dashed lines is clearly 
visible for y = 6.35. The horizontal bar indicates the region 
that was irradiated during a Ti measurement. See supple- 
mentary material for details and for other doping levels [ij- 



lattice constant jl5|. ^'^Cu NMR spectra were obtained 
by integrating averaged spin-echo signals as the frequency 
was swept, and the spin-lattice relaxation rates, T^^, 
were measured by monitoring the recovery of the nuclear 
magnetization after a saturation pulse. 

Fig. 1 shows ^•^Cu spectra for y = 6.35 and 6.45 mea- 
sured at iJ = 7 T applied along the c axis. The leftmost 
sharp line and the rightmost broad line in each panel 
are identified to arise from the Cu(l)o and Cu(l)2 sites, 
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respectively, in the CuO chains, where the subscript de- 
notes the number of the nearest neighboring oxygen ions 
along the chain [3l • While the central line at room tem- 
perature comes from the planar Cu(2) site [l6| . we find 
that it evolves in a complicated wayas T is lowered, es- 
pecially for y = 6.35 and y = 6.4 The intensity of 
Cu(2) (dotted line in Fig. 1) is strongly suppressed, and 
at the same time a new feature spO, consisting of two 
resonance lines indicated by two dashed lines emerges at 
an onset temperature Clearly, for y = 6.45, spO 

occurs at a much lower temperature than for y = 6.35. 
The rapid suppression of Cu(2) with decreasing T is at- 
tributed to the wipe-out effect, as observed in ^^Cu NQR 



in Yi_2Ca2Ba2Cu30i 



which arises from a slowing 
il, [S]. In contrast, the T- 



down of spin fluctuations 

dependence of the ^'^Cu spectra for both y = 6.5 and 6.6 
are almost identical without any signature of either spO 
or wipeout of Cu(2 ), e xcept for a moderate broadening 
with decreasing T [ij]. Interestingly, we observed the 
splitting of the Cu(l)2 line at low T for y = 6.35, which 
may suggest that spO influences Cu(l)2. In addition, it 
is noticeable that both spO and Cu(l)o broaden signifi- 
cantly at low T below 20 K for y = 6.35, suggesting the 
occurrence of glassy or incommensurate magnetic order. 
Note that such a broadening could not be measured for 
y = 6.4 and 6.45, because the magnetic order occurs deep 
in the superconducting state where the NMR spectra sig- 
nificantly weaken and become complicated. 

Fig. 2 shows the recovery of the ^'^Cu nuclear magneti- 
zation after a saturating pulse as a function of delay time 
t at various temperatures for y = 6.35 and 6.45, measured 
at the Cu(2) line. Unexpectedly, we found a step-like fea- 
ture in the data, which indicates two Ti-processes. The 
long Ti process becomes discernible at 265 K for y = 6.35 
and at 120 K for 6.45, and its fraction (i.e., 1 — ni{t) at 
the step) increases rapidly with decreasing T. Compar- 
ing with the T-dependence of the spectrum as shown in 
Fig. 1, we find that the long Tu component arises when 
the new spectrum spO becomes visible, suggesting that 
spO is due to an emerging static phase featured by the 
long Tii Since two Ti processes are evident in the 

relaxation data, we used a fitting function for magnetic 
relaxation of the central transition for 7 = 3/2 including 
two Ti components. 



m{t) = 1 - a {(1 - w) (o.le-*/'^i= + 0.9e-^*/'^i=) 
+ «;(0.1e-*/^"+0.9e-«*/^")} 



(1) 



where a is a fitting parameter which is ideally one for 
saturation recovery. Tis (Tu) are the short (long) Ti 
components, and w is the fraction of the volume governed 
by the long Tu process to the total volume. Solid curves 
in Fig. 2 are fits to Eq. (P). At low T (< 70 K), it was 
necessary to impose the stretching exponent /3 in Eq. (1) 
for the Tu recovery, which is indicative of a crossover to 
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FIG. 2: Recovery of the (normalized) nuclear magnetization 
m{t) as a function of time t for y = 6.35 (a) and y = 6.45 (b). 
In (a), m{t) deviates progressively from a single Ti-process 
with decreasing T, revealing the long Tif-process which over- 
whelms the short one at low T. Note that the short Tis- 
component is no longer detectable at 30 K. In (b), the Tu- 
process is strongly suppressed i.e., it appears below ~ 150 K 
and a large portion of the Tia-process remains down to 30 K. 
Solid lines are fits to Eq. ([T]). 



a glassy magnetic phase. 

The resulting {TuT)-\ {TuT)''^, and w as a func- 
tion of T and for all y are presented in Fig. 3. The 
results in YBCOg.e are compatible with data known so 
far [l^, revealing the spin pseudogap with the onset 
T* 145 K. The spin-lattice relaxation rate probes 
the gap solely in the spin excitation spectrum since 
(TiT)-i cc Eq^'(q)x"(q,wo) where ^(q) is the hy- 
perfine coupling constant and ujq the Larmor frequency. 
Therefore, the onset temperature of the PC and its 
doping dependence obtained by (TiT)"^ can be signif- 
icantly different from those obtained by the Knight shift, 
which probes the spin response at q = only [17|, and 
other techniques such as ARPES and optical conductivity 
which probe the charge gap [H. For YBCOe.s, T* is re- 
duced to ~ 110 K as indicated by arrows. We performed 
the measurement in another YBCOg.s crystal which has 
a longer ortho-II correlation length ^ 100 A. For this 
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sample, both p and turns out to be larger than the 
sample with shorter correlation length, being consistent 
with slightly larger T* . Upon further lowering y to 6.45, 
{Ti(T)~^ increases with decreasing T reaching a plateau 
below ^ 150 K, which is in contrast to the pseudogap 
behavior reported in a similarly doped compound [20|. 
(Ti.T)-! of YBC06.35 and YBCO6.4 is further enhanced 
with no signature of the PG as well, although Tis is not 
measurable below 100 K and 50 K, respectively, due to 
limited experimental resolution [see Fig. 2(a)]. Note 
that the T-dcpcndcncc of the data undergoes a dramatic 
change as y is reduced from 6.5 to 6.45, implying the 
existence of a critical doping just below y = 6.5. 

Fig. 3(b) shows the long relaxation rate {TiiT)~^ ver- 
sus T which arises from the spO line. For comparison, 
corresponding (TigT)"^ data of Cu(2) are also shown. 
{Tu is not detected for y = 6.5 and 6.6.) In YBCOg.ss, 
{TuT)~^ is almost constant at high T, but it starts to 
rise steeply below ~ 100 K, forming a sharp peak cen- 
tered at ~ IIK. These behaviors, together with the sig- 
nificant broadening of the spectra at low T, lead to the 
conclusion that a spin freezing or spin-glass (SG) tran- 
sition occurs at a characteristic temperature Tg ~ 11 
K for y — 6.35. We find that our data resembles 
the results of ^^Y NMR in Yi-^Ca^BaaCujOj^ 



,Sr^Cu04 21 



171 and 



i39La NQR/NMR in Laa- 
dependencies of (Ti^T)"-'^ were also observed for y 



221. Similar T- 
6.4 

and 6.45. Although we were not able to identify the lo- 
cal maximum for these doping levels due to the higher 
Tc which complicates the identification of Tg by NMR, 
one can obtain Tg ^ 5 K for y = 6.45 from the /xSR 
study of Yi-a^Caa^BaaCuaOfi which shows a very similar 
p-dependence of Tg 
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Fig. 3(c) shows the volume fraction w as a function 
of T, obtained from the fit of relaxation data m{t) using 
Eq. (1). It should be emphasized that the values of w 
themselves have no quantitative meaning, since they may 
depend on the frequency at which the relaxation rates 
were measured. Moreover, the wipeout of Cu(2) should 
lead to a significant increase of particularly at low T, 
which is indeed thought to account for the rapid increase 
of w at low temperatures [see Fig. 3(c)]. Nonetheless, the 
temperature at which it; — > 0, i.e., where the Tu process 
vanishes with increasing T should be unaffected by those 
facts. Thus, with reasonable accuracy, one can define 
the onset temperature Tg from the values extrapolated 
to w = 0, denoted by arrows. 

From these results, wc draw the p-T phase diagram in 
Fig. 4. The most striking feature is that Tq falls rapidly 
to zero at a hole concentration oip ~ 0.1 beneath the SC 
dome. At the same time, the SG transition temperature 
Tg shows also similar doping dependence, being termi- 
nated at p where Tq 0. These behaviors suggest a 
close relationship between Tq and Tg, collapsing to the 
same critical doping pc ^ 0.1. Interestingly, it turns out 
that Pc is very near the doping level at which the metal- 
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FIG. 3: (a) Nuclear spin-lattice relaxation rates divided by 
T, (TiT)"^ versus T for the short Ti^. The opening of the 
PG at T* (denoted by arrows) was detected only for y > 6.5, 
whereas there are no signatures of the PG for y < 6.45. In 
(b), (Ti^T)~^ data are presented, together with {TisT)~^ for 
comparison. In YBC06.35, (TifT)-^ forms a local maximum 
at Tg ~ 11 K, indicating a SG transition. This behavior is 
significantly suppressed for y = 6.45. (c) T-dependence of the 
volume fraction w. The onset temperature To is defined from 
the values at which — >■ (down arrows). 



insulator crossover (MIC) takes place [2J| or the Fermi 
surface is reconstructed by density- wave order [1, 0, [2^ . 
The phase diagram is also in qualitative agreement with 
that suggested by inelastic neutron scattering (INS) and 



13|,l26|. 



muon spin rotation (/xSR) studies 

A remarkable finding in our study is that the spin pseu- 
dogap is observed only near and above pc (i.e., y > 6.5) 
[see Figs. 3(a) and 4]. Note that T* of YBCOe.s po- 
sitioned near pc is much lower than the value expected 
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FIG. 4: Phase diagram of underdoped YBCOy in terms of 
hole concentration p. y values are shown on the top axis for 
convenience. AFM transition is from Ref. [l^ and dotted 
line drawn for T* is estimated from previously known NMR 
results OfiS. 



based on extrapolation of the doping dependence estab- 
lished at higher p, which is indicated by the dotted line. 
Indeed, the sudden drop of T* near pc in YBCO6..5 ap- 
pears as a crossover to the absence of T* in YBCO6.45, 
corroborating the INS results in which the suppression 
of the low energy spin excitations (i.e., spin pseudogap) 
was not detected down to 5 K in YBCO6.45 [H, [27I |28|. 
Such an abrupt suppression of the spin pseudogap below 
a critical doping level was reported in NMR studies not 
only of similarly hole-doped Yi_2CazBa2Cu30y [13] but 
also of the multi-layer cuprate Ba2Ca2Cu306(FyOi_y)2 
(29! . By the same token, the fact that pc is in the vicin- 
ity of the MIC |2J| agrees well with the disappearance 
of the quantum oscillations below ^ Pc, as well as with 
the NMR result in Bi2Sr2-rLa.rCuOfi+A that the ground 
state of the PG is metallic [30[. Therefore, we argue that 
the spin PG phenomenon may be quantum critical, stem- 
ming from the suppression of magnetism. In fact, such a 
strong competition between magnetism and the PG gives 
a good account of the observation that the PG is abruptly 
suppressed by substituting the Cu(2) sites with ^ 1% Zn 
impurities around which a local moment is induced on 
the Cu(2) sites in YBCOe.? 



31| and YBa2Cu408 [32|. 



While the nature of static order which appears at Tq is 
unclear yet, one may consider Tg as the onset of stripe- 
like charge modulation in the plane, possibly induced 
by the end Cu ion of the oxygen-filled chains, Cu(l)i. 
In fact, as effective impurities, Cu(l)i sites can cause 
Friedel-like oscillations j33l | , which may in turn induce the 
stripe-like charge modulation in the plane. This would 
be consistent with a higher Tq at lower p, where the num- 



ber of Cu(l)i is higher. The charge modulation scenario 
agrees well with signatures of broken rotational symme- 
try observed in resistivity [1] and INS [l^] measure- 
ments at p < Pc- In this case, the planar Cu(2) sites 
should be differentiated into two spatially modulated dis- 
tinct regions. In terms of stripes, the short Tig-yielding 
region is naturally related to spin stripes which consist 
of the localized Cu(2) spins, while charge stripes may 
yield the long Tu if the spin contributions to the relax- 
ation rates were almost quenched. Although the stripe- 
like charge modulation accounts for our results to large 
extent, it remains a question whether the two alternating 
regions in the stripe structure could, in practice, result in 
the (TisT)-i and {TuT)-^ values that strikingly differ 
up to three orders of magnitude. 

In summary, we performed a systematic ^"^Cu NMR 
study as a function of doping and temperature in highly 
underdoped YBa2Cu30j,, showing that static order, 
probably stripe-like, emerges at the onset temperature 
To, being followed by glassy magnetic order at Tg. The 
resulting phase diagram includes a critical hole doping 
Pc ^ 0.1, at which both Tq and Tg fall to zero. An- 
other important finding is that the spin pseudogap was 
detected only above Pc, suggesting that the spin pseudo- 
gap phase competes with the magnetic order and/or the 
static order detected at Tg and Tq, respectively. 
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